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Abstract. We discuss recent results on the cosmology of extended theories of gravity formulated in 
the Palatini approach, i.e., assuming that metric and connection are independent fields. In particular, 
| we focus on the attempts to explain the cosmic speedup with f(R) theories and on models that avoid 

the big bang singularity. The field equations for gravity Lagrangians of the form f{gp.v,R a ^fiv) 
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INTRODUCTION 

Sh . 

One of the most elementary and important lessons of Einstein's theory of gravity, 
general relativity (GR), is that the spatial sections of four dimensional space-time need 
>■ ; not be Euclidean. The Minkowskian description is just an approximation valid on 

(relatively) local portions of space-time. On larger scales, however, one must consider 
the deformations induced by the matter on the geometry, which must be dictated by 
some set of field equations. In this respect, the predictions of GR are in agreement with 
1^! ■ experiments in scales that range from millimeters to astronomical units, scales in which 

weak and strong field phenomena can be observed [1]. The theory is so successful in 
those regimes and scales that it is generally accepted that it should also work at larger 
and shorter scales, and at weaker and stronger regimes. The validity of these assump- 
tions, obviously, is not guaranteed a priori regardless of how beautiful and ellegant the 
theory might appear to the expert. For this reason, not only must we keep confronting 
the predictions of the theory with experiments and/or observations at new scales, but 
also we have to demand theoretical consistency with the other physical interactions and, 
in particular, in the quantum regime. 



Extrapolating the validity of GR at the largest scales forces us to draw a picture of 
the universe that is not yet supported by other independent observations. For instance, 
to explain the rotation curves of spiral galaxies, we must accept the existence of vast 
amounts of unseen matter surrounding those galaxies. Additionally, to explain the 
luminosity-distance relation of distant type la supernovae and some properties of the 
distribution of matter and radiation at large scales, we must accept the existence of 
yet another source of energy with repulsive gravitational properties [2]. Together those 
unseen (or dark) sources of matter and energy are found to make up to 96% of the total 
energy of the observable universe! This huge discrepancy between the gravitationally 



estimated amounts of matter and energy and the direct measurements via electromag- 
netic radiation motivates the search for alternative theories of gravity which can account 
for the large scale dynamics and structure without the need for dark matter and/or dark 
energy. 

On the other hand, the difficulties faced by GR to provide a consistent description of 
singularities and quantum phenomena at high energies (microscopic or Planck scales) 
may be an indication that we should go beyond Riemannian structures to successfully 
quantize the theory and avoid singularities. In fact, the use of the differential geometry 
of Riemannian spaces in the formulation of GR soon led to new developments that, 
in particular, put forward that there is no logical reason to require or impose at a 
purely kinematical level that the connection should be the Levi-Civita connection of the 
metric. There is no reason to assume such a strict relation between the metric and affine 
structures of the theory. In this sense, the Palatini formulation of GR was essential in 
the ADM Hamiltonian formulation of the theory and also played a relevant role in the 
development of loop quantum gravity, a non-perturbative approach to the quantization 
of GR. 

The above arguments support the view that GR may need some kind of extension 
beyond its usual formulation to address certain observational and theoretical problems. 
In this work, we review some literature related with the cosmic speedup problem and 
cosmological singularities (very early universe) framed within the so-called Palatini 
formalism, an approach in which metric and connection are regarded as independent 
geometrical objects. Though this formulation does not affect the classical dynamics of 
GR (when torsion is neglected), it does introduce interesting new features as compared 
to the usual Riemannian approach when extensions of the Einstein-Hilbert Lagrangian 
are considered. 



DYNAMICS OF PALATINI THEORIES 

We begin by deriving the field equations of Palatini theories in a very general case and 
then consider some simplifications to make contact with the literature. For a generic 
Palatini theory in which the connection appears through the Riemann tensor or contrac- 
tions of it, the action can be written as follows 




where S m is the matter action, xy represents collectively the matter fields, K 2 is a constant 
with suitable dimensions (if / = R, then K 2 = 8xG), and R a p llv = d^T"^ — d v T^ + 

r^r^g — r^P^ represents the components of the Riemann tensor, the field strength 
of the connection r"g . Note that since the connection is determined dynamically, i.e., we 
assume independence between the metric and affine structures of the theory, we cannot 
assume any a priori symmetry in its lower indices. This means that in the variation of 
the action to obtain the field equations we must bear in mind that Tg ^ T" B , i.e., we 



admit the possibility of nonvanishing torsion. It should be noted that in GR energy and 
momentum are the sources of curvature, while torsion is sourced by the spin of particles 
[3]. The fact that torsion is usually not consider in introductory courses on gravitation 
may be rooted in the educational tradition of this subject and the fact that the spin of 
particles was discovered many years after the original formulation of GR by Einstein. 
Another reason may be that the effects of torsion are very weak in general, except at 
very high densities, where the role of torsion becomes dominant and may even avoid 
the formation of singularities [4, 5]. For these reasons, and to motivate and facilitate 
the exploration of the effects of torsion in extensions of GR, our derivation of the field 
equations will be as general as possible (within reasonable limits). 
We will assume a symmetric metric tensor g^ v = g Vfl and the usual definitions for the 
Ricci scalar R = g^ v R^ v and Ricci tensor R^ v = R p ^ pv . The variation of the action (1) 
with respect to the metric and the connection can be expressed as 
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Straightforward manipulations show that dR a ^ lxv can be written as 



+ 8S m . (2) 
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where S^ v 



(T^v — r^J/2 represents the torsion tensor, the antisymmetric part of the 



connection. From now on we will use the notation P a ^^ v 
SR a p^v term in (2) in suitable form, we need to note that 
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where/' 1 = P a ^ v 8T^. Since, in general, V^y/^gJ* 1 ) = d^^/=gJ^) +2S^ 
we find that (4) can be written as 
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Using this result, (2) becomes 
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We thus find that the field equations can be written as follows 
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where pj^v] = (p a 0Mv -P a P v »)/2, T^ v = - J—^ is the energy-momentum tensor 
of the matter, and H a v ^ = — 7=^%- represents the coupling of matter to the connec- 

v s 01 v/J 

tion. For simplicity, from now on we will assume that H a v ^ = 0. Eq. (8) can be put in a 
more convenient form if the connection is decomposed into its symmetric and antisym- 
metric (torsion) parts, T" v = C" v + S" v , such that V M A V = d^A v - C® v A a - S® v A a = 
V~A V -S% v A a and V MV /=g = V£</=g -S% ay /=g. By doing this, (8) turns into 



(9) 



Example: f(R,Q) theories 

Eqs. (7) and (9) can be used to write the field equations for the metric and the 
connection for specific choices of the Lagrangian f(g^ Vl R a p^ v ). To make contact with 

the literature, we now focus on the case f(R, Q) = f{g^ v R ilv ,g^ v g a ^Ry La R v p)- For this 
family of Lagrangians, we obtain 



(10) 



where fx = dxf- Inserting this expression in (9) and tracing over a and v, we find that 
V^W^gMW] = {2^g/3)[Sl G M^ x + (3/2)s{ M^]. Using this result, the connec- 
tion equation can be put as follows 
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The symmetric and antisymmetric combinations of this equation lead, respectively, to 
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Important simplifications can be achieved considering the new variables 
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and taking the parameter a = 2/3, which implies that S x = f f , is such that S% v = 0. 
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The symmetric and antisymmetric parts of T x v are related to those of T x v by 
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Using these variables, Eqs. (12) and (13) boil down to 
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and 
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In these equations, M^ v ) = f Rg P v + 2f Q R^ v \T), and M^ v l = 2/ e 7?^ v l(r), where 
%v)(r) =%v)(f) ^R m {T)=R m {t)-\ (dpS a av -d v S a a ^. 



In the recent literature on Palatini theories, only the torsionless case has been studied 
in detail. When torsion is considered in f(R) theories, Eqs. (17) and (18) recover the 
results presented in [6]. In general, those equations put forward that when the torsion 
tensor 5^ vanishes, the symmetric and antisymmetric parts of M@ v decouple. The 
dynamics of these theories, therefore, can be studied in different levels of complexity. 
The simplest case consists on setting 5^ and R[p, v ] t0 zero > which implies that T% v = 
d v (j), where is a scalar function that can be determined by solving the equation for 
M^ v \ In GR, we find that = ln v / = g, where g represents de determinant of the 
metric. In the case of f(R, Q) theories, this result motivates the search for a rank-two 
tensor h^ v such that § = lny/^h. This allows to write the equation for M^ v in the form 
^ a (V—hh^ v ) = 0, which allows to expres as the Levi-Civita connection of h^ v . 
The relation between h^ v and g^ v can be explicitly found using (7) once the matter 
sources are specified [7]. 

The simplest solution with nonzero torsion corresponds to the case 5^ = 0. From 
(16), we see that the condition — implies that the torsion tensor satisfies — 
l(8 a v A x -8 x v A a ), where A a = S a 



R[nv] = if the connection is chosen as Tg v 



aa . With this value of the torsion, we can impose 
|A V + d v (j), where (j) is the same scalar 



function that appears in the torsionless case. It should be noted that the case S v x = 
is the most general solution of the case R\n.v] = 0- This follows from (18), which can 



be written as 



to S v aX g$ K = S^ aX g VK . Contracting this expression with 
R[fiv] 0, then the theory contains new dynamical degrees of freedom [8] 
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COSMIC SPEEDUP IN PALATINI f(R) THEORIES 

Observations of the cosmic microwave background (CMB) radiation [9], high redshift 
supernovae surveys [10], large scale structure[l 1], and baryon acoustic oscillations [12] 
suggest that the expansion history of the universe has passed through a number of phases, 
which consist on an earlier stage of rapidly accelerated expansion (known as inflation) 



followed by two periods of decelerated expansion dominated by the presence of radiation 
and dust (matter without pressure), respectively, and a current phase of accelerated 
expansion that started some five billion years ago following the era of matter domination. 
The field equations of GR in a Friedmann-Robertson- Walker (FRW) spacetime with line 
element ds 2 = —dt 2 + a 2 dx 2 filled with non-interacting perfect fluids of density p ; and 
pressure Pi, 

K K 2 d K 2 

a L 3 a 6 



+ ^ = ^P, -_ = -=r(P + *P), (19) 



where K is the spatial curvature, p = pi, and P = P, , indicate that a phase of positive 
accelerated expansion can only happen if there exists some matter/energy source that 
dominates over the others and whose equation of state satisfies Px/ px < —1/3, where 
Px and px represent the pressure and energy density of that source. A natural candidate 
to explain the current phase of cosmic acceleration is a cosmological constant A, for 
which Pa/ Pa = — 1- However, this simple proposal is hard to accept from a theoretical 
point of view. If A represents a new fundamental constant of Nature, one could expect 
new physical phenomena at cosmic scales in analogy with what happened when the 
Planck constant was discovered. If it is seen as vacuum quantum energy, then it is 
generally claimed that its observed value is too small to be in agreement with a naive 
estimation from quantum field theory, though if we apply more rigorous techniques of 
quantum field renormalization in curved space-times the predicted value turns out to be 
much smaller [13] than the observed one. For these and other reasons, there seems to be 
a widespread desire to explain the current cosmic speedup by means of some dynamical 
entity rather than by a pure constant of cosmic nature. 

When the Einstein-Hilbert Lagrangian is extended to a function of the scalar curva- 
ture, f(R) , eqs. (7) and (8) yield 
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The second equation of above implies that the connection can be solved as the Levi- 
Civita connection of an auxiliary metric h^ v = f R g^ v , while the trace of the first equa- 
tion, Rf R -If = K 2 T, implies that R = R(T) and f(R) = f(R(T)). The field equations 
for the metric g^ v can thus be written as 
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where R llv (g), R(g), and V^Vy/^ are computed in terms of the Levi-Civita connection 
of the metric g^ v , whereas R and f R must be seen as functions of T. The field equations 
of these theories in vacuum exactly boil down to those of GR with an effective cos- 
mological constant. This turned these theories into a very natural candidate to explain 



the cosmic speedup. For suitable choices of the function f(R), it could happen that the 
new gravitationally-induced matter terms that appear on the right hand side of (22) were 
negligible during earlier phases of the expansion history but became dominant at later 
times, thus allowing an expansion that closely resembles GR in the past but produces 
cosmic speedup today. One could thus explain the transition from a matter dominated 
universe to an asymptotically de Sitter accelerated one with standard sources of matter 
and radiation but without the theoretical problems posed by a strictly constant A. The 
most famous f(R) model of this kind investigated in the Palatini approach was borrowed 
from a proposal of Carroll et al. [14] in metric formalism, namely, f(R) = R — l*l 4 /R, 
where = ji 2 / K 2 represents the energy-density scale at which the effects of the mod- 
ified dynamics are relevant. Vollick [15] considered this model and showed that after 
the standard matter-dominated era, the expansion approaches a de Sitter phase expo- 
nentially fast. To see this, consider the modified Friedmann equation corresponding to a 
given f(R) Lagrangian in a universe filled with matter and radiation 
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where p m represents the energy density of the (pressureless) matter, p r is the energy 
density of radiation, and R is a function of p m only because T = — p m . In the 1/7? model, 
one finds 



.2 

R 



which recovers R K 2 p m when p^/p m <C 1 and tends to the constant value 
R vac = v^u 2 when p^j p m 3> 1. We thus see that when p^j p m <^l then (23) behaves as 
H 2 s=3 Hq R — K 2 (p m + Ap r /?))(p l i/p m ) 2 + which is virtually indistinguishable from 
GR. However, when the matter energy density, p m ~ a~ 3 , drops below the constant 

value Pji, p^i/pm > 1, then (23) goes like H 2 ps + ^K 2 p m + which tends to 

a constant and implies an asymptotically de Sitter expansion, thus confirming the late 
time cosmic speedup. 

The l/R model was soon compared with observations of type la supernovae [16], 
though such first studies were excessively optimistic about its viability. This optimism 
may have its origin in earlier studies of Palatini f(R) cosmologies which concluded 
that these theories were very poorly constrained [17], being |//y?(0)| < 10 113 one of the 
constraints coming from cosmological data. Besides the R — fl A /R theory, which repre- 
sented a small departure from GR at low matter densities, some authors also explored 
whether radical departures from the GR dynamics at cosmic scales such as f(R) = fiR n 
or f(R) = alnR could be compatible with observations. These models were confronted 
with the Hubble diagram of type la Supernovae, the data on the gas mass fraction in 
relaxed galaxy clusters [18], and baryon acoustic oscillations [19]. Though the fits to the 
data were good, the statistical analysis did not suggest any improvement with respect to 



the standard ACDM model. On the other hand, tight constraints on the family of models 
were obtained by studying the cosmic microwave background (CMB) shift pa- 
rameter and the linear evolution of inhomogeneities [20] plus the Hubble diagram of type 
la supernovae and baryon oscillations [21]. Besides finding that the f(R) = R — /J, 4 /R 
model was strongly disfavored by the data, it was found that the combined observational 
data were capable of reducing the allowed parameter space of the exponent /3 to an 
interval of order ~ 3 x 10~ 5 around /3 = 0, with a having a value similar to the cos- 
mological constant. This meant that R — aR^ R — a — a/3 \nR could be restricted to 
a tiny region around the ACDM model. More stringent constraints on this model were 
found comparing its predictions with the CMB and matter power spectra [22], pushing 
the /3 parameter to the range ~ 10 6 , thus making this model virtually indistinguish- 
able from ACDM. These conclusions have been reconfirmed by considering updated 
data [23, 24, 25] (Union and Union2 supernovae compilations plus other determinations 
[26] of the expansion rate H(z)) and strong lensing statistics [27, 28]. Causality related 
questions have also been discussed [29] in relation with this model. A different class of 
models [30], with f(R) = (R n — Rq) 1 ^, has recently been confronted with various data 
samples. The constraints on the parameters, n = 0.98 ±0.08, also place this model in the 
vicinity of the ACDM model. 

The models considered so far modify the gravitational dynamics at late times, which 
turns out to be strongly constrained by observations. Modifications at early times should 
be very weak because of the strong constraints imposed by big bang nucleosynthesis 
and CMB primary anisotropics. One could thus consider whether modifications at in- 
termediate times could be in agreement with observations. A model proposed in this 

direction [31] takes the form 1 f(R) = R + Ai#o<H^ l//(;i2Wo)2 , where H represents the 
current value of the Hubble parameter, X\ measures the magnitude of the departure from 
GR, and Xi controls the time at which the correction becomes relevant. Note that at late 
times this f(R) recovers the ACDM model (which corresponds to the limits R — > or 
X 2 — y oo). Though the background evolution of this model is not significantly different 
from the standard ACDM model for A2 = 500, 1000, which means that it can hardly be 
constrained by type la supernovae data, its effects on the CMB and matter power spectra 
are dramatic, being A2 = 1000 safely excluded. The strongest constraints are imposed 
by the matter power spectrum. This can be understood by looking at the growth equation 
for the comoving energy density fluctuations [20, 33, 34] 8 m for large momentum k 

d 2 F> k 2 r 2 

^~~^ dra ' (25) 

where* = loga(f), and c 2 = /r/ (3fR(2fRH + /r)) represents the effective sound speed 
squared. If c 2 s > 0, the perturbations oscillate instead of growing, whereas for c 1 < 
they become unstable and blow up (this happens for f(R) = R — aR^ if /3 > 0). In the 
ACDM model c 2 = 0. The form of the matter power spectrum in the exponential and 
power-law models, therefore, changes significantly with time developing an intricate 
oscillatory structure for larger k that clearly conflicts with observations, which allows to 



A similar model with a e \ R \l(^ H o) 2 Jr correction was considered in [32]. 



strongly constrain the parameter space of the models. The most optimistic constraints 
restrict the parameter A2 to the region [31] A2 > 5 x 10 4 . 

In parallel to the considerations of above, a theoretical consistency check using phase 
space analysis [35, 36] was also carried out to determine whether some families of f(R) 
models could allow for the different phases in the expansion history of the universe 
suggested by observations. It was shown that radiation, matter, and de Sitter points exist 
irrespective of the form of the function f(R) provided that the function 

does not show discontinuous or divergent behaviors. Thus models satisfying the 
condition C(R) > —3 lead to a background evolution comprising the sequence of 
radiation, matter and de-Sitter epochs. From this it follows that, unlike in metric 
formalism, theories of the type f(R) = R — j5/R n do allow for the sequence of radiation- 
dominated, matter-dominated, and de Sitter eras if n > — 1. For theories of the type 
f(R) = R + aR m — (5/R n , one finds that an early inflationary epoch is not followed by a 
standard radiation-dominated era, which conflicts with the idea that early and late time 
cosmic acceleration could be unified with this type of models [37]. In particular, for 
m > 2, the inflationary era is stable and prohibits the end of inflation; if3/2<m<3, 
then inflation ends with a transition to a matter-dominated phase, which is then fol- 
lowed by late time acceleration; for 4/3 < m < 3/2, inflation is not possible; and for 
< m < 4/3 one can have the sequence of radiation-dominated, matter-dominated, and 
late-time de Sitter without early-time inflation. 



NONSINGULAR BOUNCING COSMOLOGIES 

We have seen that cosmological observations (and local experiments too [38, 39, 40]) 
strongly constrain the form of the f(R) gravity Lagrangian at low curvatures. Though 
many f(R) models have the ability to produce late-time cosmic acceleration and fit 
well the background expansion history, they are not in quantitative agreement with the 
structure and evolution of cosmic inhomogeneities. On the other hand, it can be shown 
that the fact that matter is concentrated in discrete structures like atoms causes the 
modified dynamics to manifest also in laboratory experiments, which confirms earlier 
suspicions on the viability of such models according to their corresponding Newtonian 
and post-Newtonian properties. This is a very disturbing aspect of the models with 
infrared corrections, which demands the consideration of a microscopic description 
of the sources and prevents the use of macroscopic, averaged representations of the 
matter. A careful analysis of this point put forward the existence of non-perturbative 
effects induced by the Palatini dynamics in a number of contexts [38, 40, 41]. In this 
sense, it is worth noting that even though the ground state of Hydrogen can be studied 
using standard perturbative methods, the first and higher excited states do manifest 
non-perturbative properties [40]. Despite the fact that the modified dynamics is strongly 
suppressed in regions of high density, non-perturbative effects arise near the zeros of the 
atomic wavefunctions, where the matter density crosses the characteristic low-density 



scale of the theory and the gradients of the matter distribution become very important for 
the dynamics [see eq.(22)]. Though this certainly is an undesired property of infrared- 
corrected models, it could become a very useful tool for models with corrections at 
high curvatures. Can we construct singularity-free cosmological models that recover 
GR at low curvatures using the non-perturbative properties of Palatini theories? As we 
will see, ultraviolet-corrected Palatini models turn out to be very efficient at removing 
the big bang cosmic singularity in various situations of interest. In this section we will 
review recent efforts carried out to better understand the properties of Palatini theories 
in the early universe. 



Growing interest in the dynamics of the early-universe in Palatini theories has arisen, 
in part, from the observation that the effective equations of loop quantum cosmology 
[42] (LQC), a Hamiltonian approach to quantum gravity based on the quantization 
techniques of loop quantum gravity, could be exactly reproduced by a Palatini f(R) 
Lagrangian [43]. In LQC, non-perturbative quantum gravity effects lead to the resolution 
of the big bang singularity by a quantum bounce without introducing any new degrees of 
freedom. Though fundamentally discrete, the theory admits a continuum description in 
terms of an effective Hamiltonian that for a homogeneous and isotropic universe filled 
with a massless scalar field leads to the following modified Friedmann equation 



where p cnf 0-^Ppianck- At low densities, p / p cr i t <C 1 , the background dynamics is the 
same as in GR, whereas at densities of order p CT j t the non-linear new matter contribution 
forces the vanishing of H 2 and hence a cosmic bounce. This singularity avoidance seems 
to be a generic feature of loop-quantized universes [44]. 

Palatini f(R) theories share with LQC an interesting property: the modified dynamics 
that they generate is not the result of the existence of new dynamical degrees of freedom 
but rather it manifests itself by means of non-linear contributions produced by the matter 
sources, which contrasts with other approaches to quantum gravity and to modified 
gravity. This similarity makes it tempting to put into correspondence Eq.(27) with the 
corresponding f(R) equation 



Non-singular f(R) cosmologies 
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Taking into account the trace equation for a massless scalar, Rf R — 2f = 2k 2 p, which 
implies that p = p(R), one finds that a Palatini f(R) theory able to reproduce the LQC 



dynamics (27) must satisfy the differential equation 
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where A = ^2(Rf R - 2f)(2R c - [Rf R - If}), B = 2^R c f R (2Rf R -3f), andR c = K 2 p c . 



If one imposes the boundary condition lim/^o//? — > 1 at low curvatures, and cilqc — dp a i 
(where a represents the acceleration of the expansion factor) at p = p c , the solution to 
this equation is unique. The solution was found numerically [43], though the following 
function provides a very good approximation from the GR regime to the bouncing region 



Though the function (30) implies that the LQC Lagrangian is an infinite series, which 
is a manifestation of the non-local properties of the quantum geometry, the fact is that 
one can find non-singular cosmologies of the f(R) type with a finite number of terms. 
In fact, the quadratic Lagrangian f(R) = R + R 2 /R P does exhibit non-singular solutions 
[45] for certain equations of state [46, 47] depending on the sign of Rp. To be precise, if 
Rp > the bounce occurs for sources with w = P/p > \/3.lfRp <0, then the bouncing 
condition is satisfied by w < 1/3. This can be easily understood by having a look at the 
expression for the Hubble function in a universe filled with radiation plus a fluid with 
generic equation of state w and density p 



where Ai = -(1 +w)pd p f R /f R = (1 + w)(l - 3w)K 2 pf RR /{f R {Rf RR -f R )). Due to 



the structure of Ai, one can check that H 2 vanishes when f R — > 0. A more careful 
analysis [48] shows that f R — > is the only possible way to obtain a bounce with a 
Palatini f(R) theory that recovers GR at low curvatures if w is constant. In the case of 
f(R) = R+R 2 /Rp, it is easy to see that f R = has a solution if 1 + 2RB unce/Rp = is 
satisfied for p Bo unce > 0, where Rsounce = (1 - ?>w)K 2 p Bounce , which leads to the cases 
mentioned above. 

Besides avoiding the development of curvature singularities, bouncing cosmologies 
can solve the horizon problem [49], which makes them interesting as a substitute for 
inflation. To be regarded as a serious candidate to explain the phenomenology of the 
early universe, these theories should provide a consistent evolution of perturbations 
across the bounce, which should also be compatible with the observed nearly scale 
invariant spectrum of primordial perturbations. Investigations in this direction have 
found [50] that f{R) models that develop a bounce when the condition f R = is met 
turn out to exhibit singular behavior of inhomogeneous perturbations in a flat, dust-filled 
universe. However, since some terms in the perturbation equations blow up as f R — > 0, 
their backreaction renders the perturbative system invalid and, therefore, one cannot say 
if there is a true singularity or not. 




(30) 
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Further insight on the robustness of the bounce under perturbations was obtained [48] 
studying the properties of f(R) theories in anisotropic spacetimes of Bianchi-I type 



3 

ds 2 = -dt 2 + Y,aj(t)(dx i ) 2 . (32) 

!=1 



If one considers these space-times under the dynamics of Palatini theories with a generic 
perfect fluid, one can derive a number of useful analytical expressions. In particular, one 

finds that the expansion = and the shear a 2 = [Hi — -|) (a measure of the 
degree of anisotropy) are given by 



G 2 = P^ (Cfz + Cja+C*) 
fl 3 

where the constants Qj = —Cji set the amount and distribution of anisotropy and satisfy 
the constraint C\i + C23 +C31 = 0. In the isotropic case, Qj = 0, one has o 2 = and 
2 = 9H 2 , with H 2 given by Eq.(31). Now, since homogeneous and isotropic bouncing 
universes require the condition f R = at the bounce, a glance at (34) indicates that the 
shear diverges as ~ This shows that, regardless of how small the anisotropics 

are initially, isotropic f (R) bouncing models with a single fluid characterized by a 
constant equation of state will develop divergences when anisotropics are present. It 
is worth noting that even though a 2 diverges at /r = 0, the expansion and its time 
derivative [48] are smooth and finite functions at that point if the density and curvature 
are finite. However, one can check by direct calculation that the Kretschman scalar 
RiivapR^P = 4(Li(Hi + H 2 ) 2 + H 2 H 2 +H 2 H 2 + H 2 H 2 ) diverges at least as ~ l/f%, 
which is a clear geometrical pathology and signals the presence of a physical singularity. 
The problems when fa vanishes seem to be generic in anisotropic models of modified 
theories of gravity [51]. It should be noted, however, that the consideration of several 
fluids, fluids with varying equation of state, or fluids with anisotropic stresses [52], 
could affect the dynamics providing new bouncing mechanisms and preventing the 
extension of this conclusion to such more realistic cases. 



Nonsingular cosmologies beyond f(R) 

The previous section provides reasons to believe that Palatini f(R) models are not 
able to produce a fully satisfactory and singularity-free alternative to GR in idealized 
universes filled with a single perfect fluid with constant equation of state. Though the 
homogeneous and isotropic case greatly improves the situation with respect to GR, the 
existence of divergences when anisotropics and inhomogeneities are present spoil the 
hopes deposited on this kind of Lagrangians. To the light of these results, new Palatini 



theories were explored [48] to determine if the introduction of new elements in the 
gravitational action could avoid the problems that appear in the f(R) models. This led 
to the study of isotropic and anisotropic cosmologies of some simple generalization of 
the f(R) family in which the Lagrangian takes the form f(R, Q), with Q = R^, V R^ V , and 
assuming that R[^ v ] = 0. Using the particular Lagrangian 

f(R,R» v Rn =R + "j£ + , (35) 

where R P ~ l p 2 is the Planck curvature, it was found that completely regular bouncing 
solutions exist for both isotropic and anisotropic homogeneous cosmologies filled with 
a perfect fluid. In particular, one finds that for a < the interval <w < 1/3 is always 
included in the family of bouncing solutions, which contains the dust and radiation cases. 
For a > 0, the fluids yielding a non-singular evolution are restricted to w > jfj^' which 
implies that the radiation case w = 1/3 is always nonsingular. For a detailed discussion 
and classification of the non-singular solutions depending on the value of the parameter 
a and the equation of state w, see [48]. 

According to (7) and (9), the field equations that follow from (35) are 

/rR^iv -^gtiv + 2fQRpaR a v = K 2 r MV , (36) 

^p[V^g(fRg^ v + 2f Q Rn] = (37) 

where fa = 9r/ and /g = <9g/. The connection equation (37) can be solved in general 
introducing an auxiliary metric h a p such that (37) takes the form [y/— M^ v ] = 0, 
which implies that T p . can be written as the Levi-Civita connection of h^ v . When the 

matter sources are represented by a perfect fluid, T^ v = (p + P)u^u v + Pg^ v , one can 
show that h^v and h^ v are given by [7] 



fylV = Q (guv - ^ ~^ M M M v^ 



(38) 



h» V = ^(g^ + ^uW) (39) 



where 



a = [a 1 (a 1 -a 2 )] 1 / 2 , A = y^p7|+^ 

k±\Jk 2 -K 2 (p+P) 



Ai = v ^7gA + ^, A 2 = v ^7e 



(41) 



It is worth noting that (38) implies a disformal relation between the metrics g^ v and h^ v . 
A relation of this form between two metrics naturally arises in Bekenstein's relativistic 
theory [53] of MOND and in previous versions of it. In the MOND theory, the vector 
Uy, is an independent dynamical vector field and the functions in front of it and in front 



of g^ v depend on another dynamical scalar field. In the theory described here, on the 
contrary, the metric tensor is the only dynamical field of the gravitational sector. Note 
also that a Palatini-like version of MOND has been recently proposed by Milgrom [54]. 
In terms of h^ v and the above definitions, (36) becomes 



Rnv{h) = 
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M 



{f + 2K 2 P) j A 1? c 2 (p+P) 



2a 



Al -A 2 



Uf^Uv 



(42) 



In this expression, the functions /, Ai, and A2 are functions of the density p and pressure 
P. In particular, for our quadratic model one finds that R = K 2 (p — 3P), like in GR, and 
Q = Q(p,P) is given by 



Q_ 
2R P 



K 2 p+ f + Rp f 2\ + Rp 

2 8 JR 32 




4K 2 (p+P) 



Rp 



where f = R + aR 1 /Rp, and the minus sign in front of the square root has been chosen to 
recover the correct limit at low curvatures. In a universe filled with radiation, for which 
R = 0, the function Q boils down to [48] 



Q 



3R 2 P 
~8~ 



%K l p l_16K 2 p 



3R P 



3R P 



(44) 



This expression recovers the GR value at low curvatures, Q ps 4(fC 2 p) 2 /3 + 
32(K 2 p) 3 /9Rp + ... but reaches a maximum Q max = 3Rp/l6 at K 2 p max = 3Rp/\6, 
where the squared root of (43) vanishes. At p max the shear also takes its maximum 

allowed value, namely, a 2 ax = a/3/1 6R 3 / 2 {C\ 2 + C| 3 + C 2 l ), which is always finite, 
and the expansion vanishes producing a cosmic bounce regardless of the amount of 
anisotropy. The model (35), therefore, avoids the well-known problems of anisotropic 
universes in GR [55], where anisotropics grow faster than the energy density during the 
contraction phase leading to a singularity that can only be avoided by sources with w > 1 . 



The evolution of inhomogeneities in the quadratic model discussed here was consid- 
ered in [56], though the approximations used there to solve for the connection equation 
did not allow to see the existence of bouncing solutions. For this reason, in this case 
one cannot make any statement regarding the evolution of inhomogeneities across the 
bounce. The cosmology of f(R) and f(R^ v R^ v ) theories was also considered in some 
detail in [57]. The possibility of having a standard cosmological evolution in f(R,Q) 
models with a large cosmological constant has been considered recently [58]. 
Before concluding, we note that the f(R, Q) theories discussed here are able to repro- 
duce other aspects of the expected phenomenology of quantum gravity at the Planck 
scale [59]. In particular, without imposing any a priori phenomenological structure, the 
quadratic Palatini model (35) predicts an energy-density dependence of the metric com- 
ponents that closely matches the structure conjectured in models of Doubly (or De- 
formed) Special Relativity [60] and Rainbow Gravity [61]. This confirms that Palatini 



theories represent a new and powerful framework to address different aspects of quan- 
tum gravity phenomenology. 
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